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Abstract 



We dispel the misconception that data from SN 1987A favor the normal 
neutrino mass hierarchy over the inverted hierarchy for sin 2 #13 > 10 -4 . We 
find comparable fits for the two hierarchies. No bound can be placed on the 
mixing angle #13 even at the lex level. 
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On February 23, 1987, antineutrinos from a 2OM Type II supernova in the Large Mag- 
ellanic Cloud were detected by the Kamiokande II (KII) [l[ and Irvine Michigan Brookhaven 
(1MB) [0 detectors. The data on SN 1987A have been important in confirming the generic 
features of the core collapse model of supernovae ||. However, the average antineutrino 
energy (E Pe ) and the binding energy of the star Eb, extracted from the signals at KII and 
1MB, only marginally agree with each other [f§-|6[]. Moreover, even under the most optimistic 
conditions, the allowed values of (E Pe ) are lower than the following predictions of traditional 
supernova models 0: 

(E Pe ) = 14-17 MeV , 

(E Px ) = 24 - 27 MeV, x = //, r, 

E b = 1.5 - 4.5 x 10 53 ergs. (1) 

Including neutrino oscillations does little to resolve these discrepancies. For the favored 
large-angle solution to the solar anomaly [§], compatibility of the two data sets is not signifi- 
cantly improved [f§-fj. Furthermore, in any oscillation scenario, v e <-> v^^ oscillations result 
in a hardened initial spectrum, implying even lower (E Pe ) and increasing the discrepancy 
with theoretical predictions The degree to which the data and supernova models dis- 

agree depends on the type of neutrino mass hierarchy. There are two possibilities depending 
upon whether m 3 is the largest or smallest mass eigenstate: 

Amj 2 = m 2 — m\ > , normal hierarchy, (2) 
< , inverted hierarchy, (3) 

where oscillations between the m 3 and mi [m-i and mi) mass eigenstates are responsible for 
the atmospheric (solar) neutrino deficit. The disagreement with the theoretical predictions 
is aggravated if the hierarchy of neutrino masses is inverted and sin 2 #13 > 10~ 3 because 
the conversion of the original v e spectrum is almost complete in this case (regardless of the 
solution to the solar neutrino problem). This is a hint that the inverted mass scheme is 
disfavored by SN 1987A @||. However, it is an overstatement that SN 1987A gives a strong 
indication that the inverted hierarchy is disfavored unless sin 2 #13 < 1CT 4 ||. 

We show that the SN 1987A data do not favor one hierarchy over the other by performing 
a maximum likelihood analysis of the data from KII and 1MB in the three neutrino framework 
assuming the LMA solution to the solar neutrino problem. We perform three-parameter 
fits to T Pe , E b and sin 2 26 , 13 , where T Pe is the temperature of the electron antineutrinos 
originating in the supernova, which for a Fermi-Dirac distribution is related to (E Pe ) via 
T Pe = 180C(3)(£ Pe )/(7vr 4 ) « (E Pe )/3.15. We find that the data can be fit equally well by 
the two hierarchies with the maximum of the likelihood function in the two cases to be 
approximately the same. No constraint can be placed on sin 2 29 13 at the la level. 

Recent refinements in supernova codes suggest that the discrepancy with SN 1987A data 



may soon be resolved. The simulations of Ref. [I(J find (E Pe ) and (E Px ) below the ranges in 
Eq. ([[]). The lower T Px is attributed to the inclusion (for the first time) of nucleon-nucleon 
neutral-current bremsstrahlung which is softer than e + e~ annihilation (the other major v x 



source). A similar conclusion is drawn in Ref. [^] which explains that nucleon recoils lower 
T Px substantially. Additional energy transfer due to neutrino-matter scattering processes 
not included in supernova simulations will further soften the emergent spectra. 
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To simulate the time averaged spectrum of antineutrinos from SN 1987A, we employ a 
Fermi-Dirac distribution with equipartition of energy between the neutrino flavors. As the 
antineutrinos travel outward from their production point in the supernova, they encounter a 
density profile that falls like 1/r 3 [pT2| | . For the LMA solution, antineutrinos can experience 



at most one resonance (at density p ~ 10 3 — 10 4 g/cm 3 and characterized by Am^ and 
sin 2 2613), and only if the hierarchy is inverted. For the normal and inverted hierarchies, the 
survival probability of electron antineutrinos is given by [|TB 



P=Ple, (4) 

and 

p = P H P le + (l-P H )sm 2 9 13 , (5) 

respectively. Here, P\ e = Pi e (E u , Am^i, sin 2 29^) is the probability that an antineutrino 
reaching the earth in the v\ mass eigenstate will interact with the detector CIS db . and 
P H ~ e - sm 2 26>i3(| Amy/^a/s wm i s the hopping probability at the resonance. If sin 2 29 n 



10~ 3 , we have Ph ~ 1 and the survival probabilities for the two hierarchies are the same. 
Thus, the normal and inverted mass hierarchies are indistinguishable for sin 2 2# 13 <C 10~ 3 . 
If sin 2 2# 13 > 10~ 3 , for the inverted hierarchy p « sin 2 # 13 < 2.5 x 10~ 2 and the original 
electron antineutrinos have all been swapped for the more energetic /i and r antineutrinos 
by the time they exit the supernova envelope, resulting in a harder incident spectrum. Thus, 
the initial v e spectrum would have to be softer for the inverted hierarchy than for the normal 
hierarchy. 

In performing our statistical analysis, we fix the solar neutrino oscillation parameters at 
Am^ = 3.7 x 10 -5 eV 2 , sin 2 2# 12 = 0.79 || and the atmospheric neutrino scale |Am 2 2 | = 
3 x 10~ 3 eV 2 [15[] and fit the data taking E^, T Pe and sin 2 2^13 as free parameters. The 



time structure of the SN 1987A signal (above threshold) is shown in Fig. 1, where it has 
been assumed that the first events at the two detectors occurred simultaneously. Most of 
the events are concentrated in the first 2.75 seconds. We perform two sets of analyses to 
illustrate how sensitive any conclusion about the mass hierarchy is to the data sample chosen 
for analysis. For the first set of analyses (which we label t<i3), we use all the data reported 
by the two experiments: the 11 events observed by KII and the 8 events observed by 1MB. 
For the second set (which we call t <3 ), we only include events within the first 2.75 seconds: 
8 events at KII and 6 at 1MB. All are assumed to be antineutrino events [03]. The procedure 



employed for our likelihood analysis is that of Ref. Q up to four minor refinements listed 
below. (We refer the reader to Ref. [|J for the definition of the likelihood function and a 
description of the method). 

1. We use an improved v e — p cross section which includes Coulomb, weak magnetism, 



recoil and outer radiative corrections [17]. It is approximately 8% larger than the 
cross-section used in Ref. M. 

Instead of using the constant density approximation for the earth to calculate Pie, we 
numerically integrate the evolution equations of neutrinos through a realistic density 



profile of the earth [18] 



We approximate the time integrated spectra of neutrinos with a Fermi-Dirac distribu- 
tion instead of a Boltzmann distribution. 



3 



4. We use the detector efficiencies given in Ref. [ITII . 



To convey how these modifications (other than the different treatment of earth matter 
effects) affect the i<i3 analysis, we show the la and 2a allowed regions in (E Pe ) and T Ve for 
the case of no neutrino mixing in Fig. 0. The results from separate analyses of KII and 1MB 
data and an analysis of the combined data are shown. This facilitates a direct comparison 
with other analyses @||. The dark-shaded region is the prediction of supernova models 
from Eq. ([!]) while the light-shaded region is the result based on Ref. [ID]. The overlap with 



the predicted T Pe is significantly greater for the t<3 analysis, as can be seen in Fig. y. Also, 
KII and 1MB data from the first three seconds are consistent at the la level. 

For the inverted hierarchy, we perform three-parameter £<i3 analyses for r = T Dx /Tp e = 
1.7 (corresponding to the middle of the ranges of Eq. [I]), r = 1.4 (the lowest value of r 
from Eq. [I]) and r = 1.25 (corresponding to a v x spectrum that is softer |IO| , pT] , pO| than 



traditionally obtained), the results of which are shown in Figs. f§-|], respectively. The results 
of the equivalent t<3 analyses are shown in Figs. |7|-|9]. The figures show allowed regions at 
the la and 2a CL. in the parameter space defined by T Pe and sin 2 29^. We do not 
consider sin 2 2# 13 > 0.1 since this would violate the CHOOZ bound EI . 







E b (10 53 ergs) 


T Pe (MeV) 


sin 2 20i3 




t < 13 sec 


1.25 


3.2 


3.4 


1.3 xl0~ 6 


-41.9 




1.4 


3.4 


3.2 


1.3 x 10~ 6 


-41.6 




1.7 


4.5 


2.6 


4.0 x 10~ 6 


-41.2 




2.0 


6.3 


2.0 


1.6 x 10~ 5 


-40.6 


t < 3 sec 


1.25 


2.0 


3.6 


1.6 x 10" 5 


-35.5 




1.4 


2.0 


3.5 


1.3 x 10~ 6 


-35.4 




1.7 


2.4 


2.9 


1.3 x 10~ 5 


-35.2 




2.0 


3.0 


2.5 


3.2 x 10" 6 


-34.9 



TABLE I. Best fit values for the binding energy Eb, electron antineutrino temperature Tp e , and 
sin 2 2#i3 from three parameter fits to all the SN 1987A data and to data from the first three seconds, 
for the inverted neutrino mass hierarchy and several values of r. The corresponding logarithms of 
the likelihood function are given in the last column. 



The best-fit points are displayed in Table I, and as expected the t<3 analysis gives higher 
best-fit Tp e (since the same is true for no oscillations), and the best-fit values of sin 2 2^13 
correspond to highly non-adiabatic transitions. However, sin 2 2#i3 is not constrained even 
at the la level. Values of r greater than 1.7 lead to T Pe shifted to smaller values than for 
r = 1.7, but the sin 2 2^13 contour remains open at the la level. It can be seen qualitatively, 
that for sin 2 2#i 3 < 1CT 4 , values of T Pe from supernova codes are less inconsistent with the 
data (see the lower panels of Fig. |]-^). But, for (E Pe ) in the range found in Ref. [JTDJ] and 
r < 1.4, SN 1987A data are consistent with theoretical predictions even for sin 2 2# 13 as high 
as the CHOOZ bound. The results of Ref. [^pj indicate that both the above conditions on 
(Ep e ) and r can easily be met simultaneously. Thus, even at a qualitative level it can not 
be claimed that the inverted mass hierarchy is disfavored for sin 2 26*13 > 10~ 4 . 
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E b (10 53 ergs) 


T Ve (MeV) 




no oscillations 


3.2 


3.6 


-42.0 


t = 1.25 , P H = 


3.1 


2.9 


-42.0 


r = 1.25 , normal 


3.2 


3.4 


-41.9 


T — 1.4 , rn — U 


Q 1 
O.l 


z.o 


/ion 


r = 1.4 , normal 


3.4 


3.2 


-41.6 


r = 1.7 , Pff = 


3.2 


2.1 


-42.0 


r = 1.7 , normal 


4.2 


2.7 


-41.2 


t = 2.0 , P H = 


3.2 


1.8 


-42.0 


r = 2.0 , normal 


5.8 


2.2 


-40.6 



TABLE II. Best fit values for E b and Tp e from two-parameter fits to all the KII and 1MB 
data. Results are presented for the case in which no oscillations occur, the inverted hierarchy with 
sin 2 2#i3 = 0.01 (Ph ~ 0), and for the normal hierarchy. 

For a given value of r, the range of T Pe that is allowed for the normal hierarchy is the 
same as that for the inverted hierarchy with sin 2 2# 13 = since the survival probabilities, 
Eqs. (fj) and ([5|), are the same for both hierarchies in this limit. The T Pe for the normal 
hierarchy can be read-off from the lower panels of Figs. f§-^| at sin 2 2^13 = 10 -6 . (This 
interpretation must be made with caution, since for the normal hierarchy, it is known a 
priori that the antineutrino spectra are independent of sin 2 2#i3, which would strictly mean 
that the confidence regions should be determined for two parameters and not three). 

To make the fact that SN 1987A data did not probe the neutrino mass hierarchy even 
more transparent, we show the results of two-parameter fits in (E Pe ) and T Pe for the normal 
hierarchy and the inverted hierarchy (with sin 2 29^ = 0.01) in Table II (£<i3) and Table III 
(t<3). We do not show figures of these allowed regions since these exist in the literature |§-|6]|1 
for the t<i3 case. As has been pointed out by previous authors, if the hierarchy is inverted 
and sin 2 2#i 3 is large (i.e., Ph ~ 0), lower values of T^ e are required to fit the data than if the 
hierarchy is normal. Since in all cases the likelihoods are comparable, the data itself does 
not favor one neutrino mass hierarchy over the other. Any deductions about the hierarchy 
can only be based on the disagreement of the data with supernova model predictions, which 
are in a state of change. 

In addition to the above arguments, one must keep in mind that SN 1987A provided us 
with very limited statistics with a somewhat uncertain time sequence. Moreover, even with 
oscillations, the marginal overlap between the KII and 1MB data should suggest that any 
conclusions can only be suggestive at best. 

In summary, we have disputed the assertion that SN 1987A provides a strong indication 
that the inverted mass hierarchy is disfavored for sin 2 6*13 > 10 -4 0. The data provides no 
substantial evidence that this is the case. Rather than telling us anything about the nature of 
the neutrino mass hierarchy, SN 1987A data seem to indicate the need for more sophisticated 
supernova codes |10] which could remove the (mass-hierarchy- independent) discrepancy with 



model predictions. For a future galactic supernova, data from Super-Kamiokande and SNO 
would enable a precise determination of both Eb and T Pe [0,0] (unless sin 2 2013 > 10 -3 and 
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E b (10 53 ergs) 


T Ve (MeV) 




no oscillations 


1.8 


4.0 


-35.6 


t = 1.25 , P H = 


1.8 


3.2 


-35.6 


r = 1.25 , normal 


1.9 


3.7 


-35.5 


T — 1.4 , rjj — U 


1. 1 




-OO.D 


r = 1.4 , normal 


2.0 


3.5 


-35.4 


r = 1.7 , Pff = 


1.7 


2.4 


-35.6 


r = 1.7 , normal 


2.4 


3.0 


-35.2 


t = 2.0 , Pff = 


1.8 


2.0 


-35.6 


r = 2.0 , normal 


3.1 


2.5 


-34.9 



TABLE III. Same as Table II but only data from the first 2.75 seconds are used. 



the hierarchy is inverted f22|) and determine the sign of Am^ if sin 2 2#i3 > 10 4 ||l"3| , |22 
independently of supernova models. 
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FIG. 1. The time structure of the SN 1987A signal. 14 of the 19 events occurred in the first 
2.75 seconds. 



8 



53 15 
o 



KII 



7TT r 
1 1 

'<< 

r< 

i ; \ imb 

i i \ \ 
i < \ \ 

\ i i \ \ 
\ i i \ \ 
\ \ i i \ \ 

1 \ M I \ \ 



\ - i 1 x \ 




T- (MeV) 

FIG. 2. The la and 2<t allowed regions for and T^ e from a two-parameter fit to SN 1987A 
data in the case that no oscillations occur. The results from separate KII and IMB analyses 
(dashed), and a combined analysis (solid) of the data sets are shown. All data in the first 13 seconds 
are included. The dark-shaded region is the range of values from Eq. (fj) and the light-shaded region 
is based on Ref. . 
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FIG. 3. Same as Fig. but only data from the first 2.75 seconds are included. 
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FIG. 4. la and 2a allowed regions from a three-parameter i<i3 analysis of combined KII and 
1MB data for the inverted hierarchy. A ratio r = T Dx /T Ue = 1.7 is assumed. 
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FIG. 5. Same as Fig. | but with r = 1.4. 
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FIG. 6. Same as Fig. | but with r = 1.25. 
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FIG. 7. la and 2a allowed regions from a three-parameter t<3 analysis of combined KII and 
1MB data for the inverted hierarchy. The ratio r = Tp x /T Ue is 1.7. 
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FIG. 8. Same as Fig. @ but with r = 1.4. 
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FIG. 9. Same as Fig. @ but with r = 1.25. 
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